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STUDIES IN THE PHYSICOMATHEMATICAL THEORY 
OF ORGANIC FORM 


N. RASHEVSKY 
THE UNIVERSITY OF CHICAGO 


In continuation of previous studies the theory is developed on the 
assumption that the form of any organism is determined by require- 
ments to perform definite biological functions. A previously outlined 
theory of the form of plants is developed further, showing how the con- 
ditions of mechanical strength together with the specifications of the 
total mass and metabolism, may quantitatively determine not only the 
general form of the plant, but even the number, size, shape and shades 
of the leaves. Next the form of animals, as required by mechanical 
conditions and bythe different types of possible locomotions is discussed. 
A mathematical theory of locomotion of snakes in relation to their shape 
is outlined. Next is discussed the form and locomotion of quadrupeds. - 
A number of theoretical relations, which describe the shape of an ani- 
mal, are derived and compared to available observations. After that the 
theory of flight of birds and insects is discussed, and again some form 
relations comparable with observations, are discussed. Finally a set of 
equations is outlined, which determines not only the external shape, but 
also the internal structure of animals. Different relations pertaining’: 
to some ‘inner organs are derived and compared with available observa- 
tions. The paper ends with a brief discussion on the shape of unicellular 

- organisms. 


Introduction 


The problem of organic form is justly considered as one of the 
most difficult problems of biology. Rather pessimistic views as to the 
possibility of a physicomathematical treatment of that problem have 
been expressed not only by some philosophically minded biologists 
(L. v. Bertalanffy; 1928), but even by such physicochemically ori- 
ented investigators as A. Meyerhoff (1926). A definite doubt as to 
whether the problem of form lies at all within the scope of physics 
is felt through many of the arguments. The physicist of course knows 
that such skepticism is not justified. Not only does crystallography 
offer a problem of form in physics, but whenever a physicist deals 
with fields, such as diffusion fields, electromagnetic fields or others, 
questions of spatial relations, and hence of form, enter explicitly. © 

In view of the importance of the diffusion fields in biology, N. 
Rashevsky (1929) has indicated a mathematical approach to the prob- 
lem of organic form by considering the combined effects of diffusion 
field and surface tensions. Later on the theory was amplified by the 
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introduction of diffusion drag forces (Rashevsky, 1938, 1940) which 
threw some light on the form of single cells. 

A rough attempt to develop a physico-mathematical theory of the 
form of metazoan animals was made by Rashevsky (1933), who ar- 
rived at a gross theoretical classification of some existing forms. Sub- 
sequently he proposed an improved theory (Rashevsky, 1938, Chapter 
xiv), which so far remained only an outline. The principal reason for 
that was a realization that although a reduction of the problem of 
organic form to physico-chemical interactions may be possible in 
principle, the working out of actual specific instances would be so 
complex, as to make the task technically hopeless. Therefore recently 
Rashevsky (1948a) proposed.a different approach to the whole prob- 
lem, seeking to solve it by the introduction of one or a few simple for- 
mal principles or postulates, which need not immediately be reduced 
to a physicochemical mechanism. Such a procedure is not unfamiliar 
to a physicist. 

The postulate proposed by Rashevsky is that the geometric form 
of any organism is determined by the requirements of certain bio- 
logical functions. Since these are almost innumerable, it is proposed, 
for the beginning, to limit the theory to a few of the simplest ones. 
As such, the metabolic level, the speed of locomotion and the amount 
of nervous coordination have been suggested. Given an organism of a 
mass M , average metabolism q ergs per unit mass, average amount 
of nervous tissue nM , the average velocity of locomotion v will be 
determined by the geometrical structure of the organism. In general 
however, the same results may be obtained by different structures. 
Therefore, to make the problem definite, Rashevsky introduces a new 
principle, namely, that of all the structures which give us the desired 
functions, the simplest one actually occurs in nature. This of course 
requires a proper definition of simplicity. 

In the case of plants, in which usually we may put v = 0 and 
n = 0, the form is to be determined by the considerations of mechani- 
cal strength and of sufficient size of food conductors, which enable 
the plant to have a given metabolic rate q for a given mass M. . 

In the above mentioned papers, Rashevsky has outlined very 
roughly how the theory is to be developed. In case of animals, terres- - 
trial locomotion without extremities is limited as to speed. Only ani- 
mals with low velocity of locomotion, will be possible without extrem- 
ities (snakes, worms). In that case the ratio of length to width of 
the animal and hence its shape is determined by M »qand v. High- 
er velocities of locomotion may be obtained with either four or two 
extremities. Though the second case is structurally simpler, it re- 
quires much more nervous coordination, since standing on two legs. 
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is mechanically an unstable process. Hence, where the necessary co- 
ordination is not sufficient for maintaining a posture on two extrem- 
ities, quadruped forms result. The size of extremities is determined 
by M,qandv. 

In this paper we propose to elaborate further the above sugges- 
tions and derive more relations, which determine the shapes of plants 
and animals. It must however be emphasized, that in this study we 
still remain largely on a purely abstract, theoretical level. The actual 
situations are much more complex. As in other branches of mathe- 
matical biophysics, we must start with idealized, oversimplified cases 
and study at first all conceivable situations though some of them may 
never occur in nature. Only when a great many of such conceivable 
situations have been studied will a comparison of the theory with 
observations indicate which of the situations actually occur. 


Although in deriving any relations from too abstract and over- 
simplified assumptions we should not expect those relations to hold 
actually, nevertheless, even in such cases a mathematically derived | 
relation may at least very roughly approximate reality. There- 
fore, wherever possible, we shall make comparisons between theory 
and observations. We must however constantly keep in mind, that 
whenever such comparisons show any agreement between theory and 
observations, this is not to be considered as constituting any proof 
of the theory, but rather as an illustration of what the theory may do 
and as a stimulus for further theoretical and experimental research. 

It is of course impossible to cover at this stage the whole field 
of organic form. We shall therefore confine ourselves to a few se- 
lected topics. In section I we shall give a rather abstract theory 
of the form of plants, extending the previous studies (Rashevsky, 
1943d). Section II deals with some types of locomotion of snakes in 
relation to their form. In section III we pass to the consideration of 
quadrupeds. From considerations of mechanical strength we shall 
derive two relations which partly determine the form of the animal. 
In section IV, V and VI, the terrestrial locomotion of quadrupeds and 
bipeds is discussed in relation to their form, and another relation, 
determining the shape of the animal is derived. Section VII deals 
with the flight of birds and insects as related to their general shape. 
Finally in section VIII the relations of some physiological functions 
of the organism to its internal structure are discussed and several 
quantitative expressions derived, which determine that internal struc- 
ture. Thus not only the external form but also the internal structure 
of animals is taken within the scope of this study. We conclude with 
a brief discussion of forms of unicellular organisms. 
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I. Forms of Plants 

It has been shown by N. Rashevsky (1943a) that to some extent 
the general shape of a plant may be determined by the mass M and 
metabolism gM, principally from mechanical considerations. Denot- 
ing by I) the length of the main trunk, by 7) its radius, by / and r the 
average length and radii of the branches, by n the total number of 
branches and by 8 the average density of the plant, Rashevsky estab- 
lished six general equations for the determination of those six quan- 
tities in terms of gq and M. By making specific assumptions about 
some functions, which were left originally indeterminate, Rashevsky 
(1943d) arrives at explicit expressions for the above six quantities. 
But five of those largely determine the gross shape of the plant. Small 
l, and 7, and large n, give us a bush. Small 1, , large 7), small n and 
large 1 correspond to a form of an apple tree, while large J, , 7), and 
n but small l and r produce a tree like a poplar (Rashevsky, 1943d). 
Two pathways are open for further studies. One is to refine the 
specific assumptions made by Rashevsky (1943d) and compare the 
theoretical results with possible measurements on plants. The other 
is to remain on the abstract level, and to investigate to what extent 
further details of the structure of a plant may be determined by M 
and q alone. We shall here follow the second pathway. 


It is very doubtful whether actually M and q are the only inde- 
pendent variables, which determine the shape of a plant. Actually, 
the metabolic processes differ in different plants, and q represents 
therefore a sort of an average quantity. More detailed specifications 
of the qualitative character of the total metabolism will actually have 
to be given. Environmental factors also probably play a role. Since 
we remain at present on a purely abstract theoretical level, it is of 
interest to show to what extent M and q alone may determine the de- 
tails of a plant’s external structure. 

Besides the variables introduced above, let us introduce the fol- 
lowing additional ones: A,—the average length of a leaf, A.—its aver- 
age width, m—the total number of leaves, d—the length of the leaf- 
stalk, 7,—its radius, r,—the average radius of the vein of a leaf, and 
n'—the average number of veins per leaf. For simplicity let us con- 
sider the thickness of a leaf as a constant, so that its mass is approxi- 
mately proportional to 4,4,. We also shall at first consider only 
branches of first order. Higher order branches will be introduced 
. subsequently. 


As before, (Rashevsky, 1943d) we have the following equation 
for the total mass: 


M = 28 (lpr? + nlr?)., (1) 
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The total metabolism being assumed to be proportional to the total 
area of the leaves, we now have, denoting by k a constant 


qM =kA,A.m. (2) 


The length of a branch is limited by its mechanical strength and its 
mass. The latter is a function of the size of the branch and of the 


mass of all the leaves on the branch. We therefore have a relation of 
the form: 


— m 
=F(r,8,7- 4,42). (3) 


Mechanical considerations give us again the relation 
lL =fo(,5,M). (4) 


The requirement that the cross-section of the trunk and of each 
branch should be sufficient to allow ample flow of metabolites gives 
(Rashevsky, 1943a, d): 

ro =f1(qM , 8), (5) 


M 
r=h(—, 3). (6) 


The number of leaves along the axis of a branch is limited by the 
width A, of the leaf and is of the order of //4,. The number of 
leaves that may be arranged along the circumference of a branch is 


A 
also limited by 4.. It is of the order of 2a(r + d + As, since 


2n(r + d + A,/2) is the length of the circumference along which the 
centers of the leaves are arranged. Hence denoting by a, a coefficient 
of proportionality, the total number m/n of leaves per branch is of 
the form 


m/n = 2na,(r + d+ A,/2)1/A2?. (7) 


The distance between the branches along the trunk is of the or- 
der of magnitude of d + A,, or is proportional to it. The total num- 
ber of branches that can be arranged on the same level of the trunk 
is approximately 27,/r. Hence the total number of branches is ap- 
proximately 


R—arolo/t(d + A,).; (8) 

A consideration similar to the one which leads to equations (5) 

and (6) gives: 
i ahs %, = fo(qM/m) ; - (9) 
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while consideration of mechanical strength leads to a relation: 
G—faltg seis) s (10) 


The cross section of each vein must be such as to permit the flow of 
necessary metabolites. The metabolism per vein being gM/mn' we 
have a relation 


r= f.(qM/mn). (11) 


Since the metabolites flow into the leaf through the surface of the 
vein, which is approximately 2ar,4., since A./2 is approximately the 
length of a vein, that surface must be of such a size as to permit the 
proper flow. The principle of maximum simplicity suggests that that 
surface be a constant multiple of the size necessary to permit the 
proper flow. We therefore have a relation: 


M 
ee = fx (Ae Ts). (12) 

mn 
Finally the mass of all the veins may be assumed not to exceed 
a certain fraction of the total mass of the leaf. Or we may introduce 
again considerations of mechanical strength, since the veins form 
part of the supporting framework of a leaf. In either case we get a 

relation of the form 


fe(A1, 42,73,”') =0. (13) 


Altogether we have (13) equations for the determination of the 
13 varieties: lei; 7o503t5 by iy © is Zia He he te ee 

If we assume that the metabolism of a leaf is proportional to its 
chlorophyll content, then we must consider the constant k in equation 
(2) as proportional to that content. We may leave k as a third inde- 
pendent variable, which by the way, determines the lighter or darker 
shade of a leaf. With proper choice of the different functions involved 
in the above equations, it may happen that those equations will have 
solutions only for values of k , lying within a relatively narrow range 
k, <k <k,. This may account for the circumstance that we do not 
have a too wide variation of shades in leaves. Saprophyte plants will 
have to be treated separately, and do not fall within the scope of this 
theory. 

It has been suggested (Rashevsky, 1943a) that the metabolic 
rate of a leaf may determine its finer shape. If for a large k, and 
large A, and A;, n' turns out to be small, the leaf may not receive 
enough metabolites at the remotest points. Therefore, while for large 
values of 7’ a form of Figure 1a will occur, for small n’ only the form 
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of Figure 1b will be possible. Hence large k, A, and A,, but small n’ 
give us dented leaves, while small k, A, and A,, and large n’ give us 
smooth-edged leaves. But with a large k and very large n’ the leaf 
will have a smooth edge, while small & and very small n’ wiil give a 
dented edge. 


| 
a b 
FIGURE 1 


As emphasized above, the assumptions underlying this outline do 
not need to be true. The important result is that a variety of plant 
forms may conceivably be determined in great detail by considera- 
tions of relatively simple biological functions. The next step to take 
is to make specific assumptions about different functions which have 
been left indefinite, and to see to what extent explicit expressions 
thus obtained will agree with observations. After that we should in- 
vestigate theoretically the effect of the individual metabolic reactions, 
thus introducing instead of one q, several independent variables. In- 
stead of a single average 5, several such quantities may have to be 
introduced for the different parts of a plant. Thus step by step we 
may approach a theory that will cover a number of actual observa- 
tions. 

We shall conclude this section by indicating how branches of 
higher order may be taken into account. For simplicity let us limit 
ourselves to considerations of branches only, as has been done before 
(Rashevsky, 1943a, d), omitting the problem of leaves. The above 
considerations which comprise the role of the leaves may then be 
readily combined with the following ones. 

Consider n, first order branches and 7, second order branches. 
Let l,, b , 7:, 72 be their corresponding lengths and radii. 

We now have: 


M =28 (br? + MhL72? + Nlre?), (14) 


and 
gM =k(mbr, + Nelet2). (15) 


Each first order branch carries n,./n, second order branches. 
Hence the weight of each first order branch plus all second order 
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branches which it carries, is 
Ne 
ns (7? ain =" loo?) ° 
Ny 
For mechanical strength we must have: 
Ne 
r,=f[8, (L772? + — 17,7) ] ’ (16) 
1 


and similarly 
1, = f (8,72). (17) 


We also have equations corresponding to equations (4) and (5). 
Moreover we have 


qM 


1, = fil , 9) (18) 


‘1 
and 


M 
n=fl— 


2 


, 8). (19) 


Altogether we have eight equations for the nine variables: §, 
lo, To, M1, 1%1,h., M2, 12,1. The ninth equation is provided by the re- 
quirement that the length 1, of the second order branches should be 
of the order of the distance between the first order branches (Figure 
2). The number of branches of first order that can be placed at the 


FIGURE 2 


same level around the circumference of the trunk is 
2a) / 274 = R/T: . 
If x is the distance along the trunk between first order branches, then 
we have: 
My lo 
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putting 1, = a.x , a» being a coefficient of proportionality, we find: 


ANP oly 
ly — ? (20 ) 


which gives us the ninth equation. 

If we consider third order branches, we introduce three new 
variables n; , 7; , ls, the number, length and radius of the third order 
branches. This introduces additional terms into equations (14), (15) 
and (16). Moreover, we obtain two new relations for 7; , which cor- 
respond to equations (17) and (19). And a third new relation is ob- 
tained by postulating that the J, should be a given fraction of the dis- 
tance between second order branches. Thus we shall have sufficient 
additional relations for the determination of the additional variables. 
The process may be extended to branches of any order. 


II, Locomotion and Form of Snakes 

There are apparently at least three different types of locomotions 
of snakes (Mosauer, 1931, 1932). Of those we shall consider the one 
which is regarded as the most common one. According to Mosauer 
the snake assumes an approximately sinusoidal shape and moves for- 
ward in such a way, that each point along the body follows the same 
sinusoidal path. The physical mechanism of such a locomotion seems 
at first somewhat puzzling. As we shall presently see, however, this 
type of locomotion may be obtained by assuming that the coefficients 
of static friction of the snake against the ground are different in the 
direction of the body axis and in the direction normal to it. The rea- 
sons for such a possibility will be discussed presently. 

Inasmuch as the exact shape of the curve formed by the body 
of the snake has never been ascertained, and may in a strict mathe- 
matical sense be not a sinusoid at all, we shall in the following dis- 
cussion make some very crude approximations, by assuming instead 
of a smooth wavy line a zigzag line (Figure 3), of total length L. The 
expressions derived are more of the nature of dimensional relations 


FIGURE 3 
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and establish orders of magnitudes, rather than exact equations. 
However, even those rough approximations throw interesting light on 
the problem. Let the length of each straight segment be 2. Let at 
every point a muscular force f per unit length tend to displace the 
different parts of the body at a right angle to AA’, as shown by the 
arrows (Figure 3). The component of this force along the snake’s 
axis iS 


fi=fsin ¢; (1) 
the component normal to the direction of the body axis is 
fr—feos¢. (2). 


Let the static friction 7; for displacement along the axis be less than 
the static friction 7, for displacement perpendicular to that axis: 


Pieerp e (3) 


This may be due to either a peculiarity of the surface of the skin 
of the snake, or to the fact that the snake, due to its weight, makes 
an impression in the ground. Such an impression tends to prevent 
lateral displacements. Longitudinal displacements may be easier, be- 
cause of the wedge-like action of the snake’s head. The total force 
acting along the snake’s axis is 


F,=f, L=fLsin¢. (4) 
The total force acting normally to the axis is 
F,=jnl —={L cos ¢. (5) 
If 
Eat | Jeera ge, (6) 


then the body of the snake will glide along its own axis. Let the ve- 
locity of this gliding be v em sec". 

A zigzag, or a wave-line shape of the body of the snake is pro- 
duced by contraction of muscles on the “concave” sides of each 
“peak”, and by relaxation of the muscles on the “convex” side of it. 
Hence if the body of the snake glides along its axis with the velocity 
v, a wave of contraction must move along the body of the snake with 
the same velocity v. If the zigzag line is periodical, as we shall con- 


sider to be the case, then, denoting by 24 the wave length, this con- 
traction wave will have a period 


P=27/0. (7) 
or a frequency 


y=v/2). (8) 
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The dry friction of the moving snake against the ground being 
independent of the velocity, any force acting along the axis, which is 
strong enough to set the body in motion, will keep this motion accel- 
erated indefinitely, if other limiting factors do not enter into play. In 
the present case such a limiting factor is found in the condition, that 
as the snake moves along its axis with the velocity v, a contraction 
Wave must move along its muscles with the same velocity. A snake 
can not move with respect to the ground with a velocity larger than 
the velocity of the contraction wave with respect to its body. Thus 
the factors limiting the velocity of the snake must be sought in the 
properties of the snake’s muscular apparatus. 

From (8) we have: 

w= 2iv. (9) 


Since 2A < L, the velocity v of the contraction wave is limited fun- 
damentally by ». This latter is probably limited by the viscosity of 
the muscle (A. V. Hill, 1927). The greater »y , the faster must be the 
velocity u of contraction of the muscle, and that velocity is limited 
by the viscosity. 

As an illustration we may consider the following line of thought. 
For a viscous rod of length /, cross section s?, viscosity 7 , we have 


1dl 1F 
ee (10) 
Ldt 3n 8s? 
where F is the total force. Hence, if a denotes a constant and 
i= as*, (11) 
as happens to be approximately the case with muscles, then 
Clee at (12) 
1=—-=>-—. 
dt 3y 


The greater the amplitude A of the wave, the larger must be the rela- 
tive shortening Al of the muscles. It is also larger, the greater the 
thickness s of the body of the snake. Moreover, the smaller the wave- 
length the larger the curvature of the body and therefore the larger 
Al. Hence 


Al «<As/i, (13) 
and 
Al As 
Uun—_—-=—rvAl « —y». . (14) 
T A 


Since the length J of the contracting mass of muscles is of the 
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order of 1, 
lw,, (15) 
hence, from (12): 
U x ad/3n . (16) 
Equations (14) and (16) give 
Asy « af?/3y , (17) 
or 
y « ad?/8Ans . (18) 
Introducing this into (9), and making the plausible assumption that 
A ods; (19) 
we have 
v « 2aL'/3Ans . (20) 


The force f per unit length may be assumed to be proportional 
to s?, because of its approximate proportionality to the cross section 
of the muscles. If we increase s , we increase both the circumference 
and the thickness of the muscles. Hence, their cross section varies 
approximately as s?. Hence, denoting by a, a constant 


f=a,s? (21) 
Therefore (4) and (6) give 
a?Lsing271, (22) 
or 
. = Ti T) 
sing =— 23 
? 3b. @,V so) 


where V is the volume of the snake. 
But, from Figure 3, 
Awilsing, (24) 
hence 
Aw {nr /a,V. (25) 
Introducing (25) into (20) we find, remembering (15) and 
(19): 
vw L?V/nsr, . (26) 
The quantity v is the velocity of the snake along its axis. Along AA’ 
the velocity is 
0=VCOS¢. (27) 


aa the density of the snake is close enough to one, we have from 


N. RASHEVSKY 13 


: rT 
sin 6 =—.,, 
oar (28) 
hence 
T;? 
cos¢= {1 — : (29) 
a,?M? 
Therefore (26) and (27) give 
= 0M. fa 
Vea ahs ; (30) 
NST; a,2M? 
Since 
L?=M, (31) 


s may be eliminated from (30), and this determines L in terms of 
M and V; ss is then determined from (31). Hence v and M, deter- 
mine the shape of the snake. 

It seems paradoxical that v should increase with M. If ™ in- 
creases itself with M , the latter may drop out of the equation. 

The weakest points in the reasoning are relations (13) and (15). 
In (13) we may perhaps consider instead of 4 the length of the indi- 
vidual muscles. This will introduce a parameter characterizing the 
inner anatomical structure. Instead of (15) a more exact relation 
should be established. 

A more rigorous procedure would be to establish from more ex- 
act considerations of the viscous losses a maximum value of » in 
terms of 1, a relation to replace (18). Let the relation be 


y—v(A) 
then the exact equation (9) gives 
y=2iy(d). (32) 


It may happen, that 4y(A) has a maximum for 4,, << L. In this 


case the equation 
dfdy(dA)]/di=0 | (33) 


may be considered as determining the actual value of 4. 

The physical reason why v increases with M is due to the fact 
that with the assumption made, the total force is proportional to M 
[equations (4) and (21)], while the viscous resistance increases with 
L , and decreases with s. 

It must also be kept in mind, that the above considerations hold 
only for uniform motion. The acceleration at the beginning will be 
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expressed in a more complicated ‘way, and it seems likely that, at 
least approximately, the acceleration will not depend on M. 

If we wish to connect the velocity of locomotion to the available 
metabolic energy Mq, we must compute for a given » and /, and 
hence v, the energy losses E per unit time, and also remember that 
the energy loss due to friction against ground is 7v erg sec”. Hence 


E+ 7r,v=Mgq. (34) 


If E is known as a function of v, equation (34) determines the 
latter in terms of M and q. 

Since the actual values of 7; and 7, , as well as of other constants 
are not known, it is impossible at present to compute the limiting 
speed of locomotion for a given q. It appears plausible, however, that 
due to the very large frictional area, those speeds would be low. Ac- 
cording to Mosauer (1931) the largest observed speeds of snakes 
are of the order of 1 meter per second. To obtain higher speeds, dif- 
ferent bodily structures must be used, involving extremities. 

A further improvement of the above outlined theory would con- 
sist in considering instead of a zigzag line a continuous line. One 
might then investigate what kind of continuous periodical line gives 
a maximum velocity, other parameters being kept constant. In any 
case the possible relation between the size and shape of snakes and 
their speed should be made the object of careful quantitative studies. 


III. Some General Considerations on the Shape of Quadrupeds 

In a very rough way, a quadruped may be considered from a 
mechanical point of view as a bar supported at its ends. This sets 
definite limits to the length of the trunk for a given width. The 
actual supporting structure is the vertebral column. However, not 
only the ribs, but even the soft muscular and other tissues contribute 
to some extent to the mechanical strength of the trunk. The study of 
different factors and conditions affecting the mechanical strength of 
the trunk of quadrupeds is a most interesting, but extremely complex 
problem. Its development will result in many relations which will 
offer clues for biometric studies. 

We shall not tackle that general problem at present, Inasmuch 
as we are as yet confining our study to extremely oversimplified pic- 
tures, we shall, for the sake of illustration only, investigate the situa- 
tion that would prevail if the trunk of a quadruped could be consid- 
ered as a uniform bar of length L and width W , supported at its ends. 

If such a bar is uniformly loaded with a load p per unit length, 
then the sagging f in the middle is (Geckeler, 1928, p. 166) equ. 16): 
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f < pLt/W?, (1) 
hence the relative sagging 
f/ Lop le Ww. (2) 
By definition 
p=M/L, (3) 


where M is the mass of the trunk, which approximately may be set 
equal to the total mass of the animal. Hence 


f/L « ML?/W+. (4) 


The relative sagging should not exceed a certain threshold value, 
lest the bar breaks. Assuming this value to be a constant, and that 
7/L is always a given fraction of that constant, we have: 


ML?/W+ = Const. or L?/W4 « M". (5) 
Moreover, we have approximately 
LW7=M. (6) 


Expressions (5) and (6) give 


L 
16, x Wee Ww foe ene ae M-'8; W o Lee (7) 
The heavier the animal, the shorter its trunk relative to its width. 


fete) 
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Ficure 4. Width W of trunk plotted against length L of trunk on logarith- 
mic scale. The straight line has the theoretically predicted slope 3/2. 1 — cat; 
2 — dog; 3 — dog; 4 — goat; 5 — ass; 6 — saw; 7 — elk; 8 — ox. Data from 


Muybridge, (1887). 
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Of course we cannot expect such simple relations to hold exactly 
in reality. It is interesting, however, that to some degree of approxi- 
mation the relation W « L*/ does actually hold, as seen from Figure 
4, where W is plotted against L on logarithmic scale. We must how- 
ever remember what has been said in the Introduction in regard to 
such agreements. 

Similar crude considerations lead us to expect that the cross Sec- 
tion of the legs will be roughly proportional to the mass M of the 
quadrupeds, since they have to be mechanically strong enough to sup- 
port the body. Such an approximation is even cruder than the one 
made above about the trunk. The strength of the legs is determined 
not only by their role as “pillars” supporting the body, but by numer- 
ous other functions. The trend towards thicker legs in heavier ani- 
mals is however unmistakable. Actual data are very scarce. To make 
even a rough comparison, the following procedure was used. Denot- 
ing the thickness of the extremities by d, its cross section is roughly 
d?. We expect approximately d? « M, or because of relations (7), 
d? « L¢+ ord « L?. Hence d/L « Lor 


ad/L « MM, (8) 


The values of d/L were taken from pictures in volume 2 of the Anat- 
omy of Vertebrates by R. Owen (1866) ; the values of M from tables 
compiled by G. Crile (1941, p. 284) and from C. Oppenheimer and L. 
Pincussen (1930, p. 732). The thickness of the femur was taken for 
d, the length between the shoulder blades and the pelvis for L. The’ 
vaules of d are extremely inaccurate. The results are plotted on Fig- 
ure 5. 


100 200 300 500 1000 
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FIcuRE 5. Ratio of thickness d of extremities to length L of trunk, plotted 
on logarithmic scale against mass M. The straight line has the theoretically 
predicted slope 1/4. 1 — lion; 2 — polar bear; 3 — horse; 4 — ox; 5 — giraffe; 
6 — hippopotamus. For sources of data see text. 
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Thus from a purely theoretical point of view we see that the 
mass M of a quadruped in principle may determine the size and shape 
of its trunk and the thickness of the extremities. This already speci- 
fies to a large extent its general form. 

In the next sections we shall see that an additional specification, 
namely the length of the extremities, may be obtained from considera- 
tions of locomotion. 

It must be remarked that relations (7) determine L and W , and 
hence the surface of the animal in terms of M. Since Rubner’s law 
requires that the total metabolism varies approximately as the sur- 
face of the body, this raises the question as to whether M and q may 
be considered as independent. However the surface of the extremities 
is not negligible as compared with that of the trunk, and therefore 
no definite conclusion can be drawn at this stage. Here is however a 
problem to be kept in mind. 


IV. General Theory of Quadruped Locomotion 


This problem shall also be treated in a very rough way. The 
idealized quadruped shall be assumed to be a bar of length L, sup- 
ported on four extremities, all of the same length /. In order of de- 
creasing simplicity we shall first discuss the case, when the extrem- 
ities have no other joints, except those connecting them to the trunk. 
Subsequently we shall consider extremities with other joints. We 
shall neglect the lateral width of the trunk and the small lateral move- 
ment of the whole body, which results during the locomotion. Things 
will be considered to happen so, as if the trunk and extremities were 
infinitely thin. We shall also consider the simplest possible time se- 
quence of the movements of the extremities, namely a rear and a 
front extremity move simultaneously forward or backward, remain- 
ing parallel. Whether the pair of simultaneously moving extremities 
is contralateral or collateral, is immaterial for our considerations be- 
cause of the assumed infinite thinness of the body. The “idealized 
quadruped” has thus at any moment the shape represented on Figure 
6. 

The theory of any motion of such a system forms a part of the 
general dynamics of chains of linked levers. The general equations 
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for such system have been developed by O. Fischer (1906). [C.f. also 
Winkelman and Grammel, 1927, p. 453]. Unfortunately, the non- 
linear differential equations involved are practically insoluble even 
in the simplest cases. While Fischer’s theory in its general form, has 
been applied to some simple cases of systems of two levers. only, it is 
at present of more interest from a mathematical point of view rather 
than from a biological. Some very drastic approximation method _ 
must be developed before the mathematical theory can be success- 
fully applied to complex practical problems. Subsequently we shall 
suggest such an approximation method. For the present we must re- 
frain from the mathematical rigour of the already developed theory 
of linked levers, and attempt to obtain specific, even though crudely 
approximate, results. The further development of the theory of ani- 
mal locomotion depends however entirely on some practical approxi- 
mation method to be adapted to O. Fischer’s theory. 

We shall make here a remark on the relation between biped and 
quadruped locomotion with rigid jointless extremities. It is readily 
seen that a biped could not have any locomotion at all in such a case. 
In order to impart a movement to the extremity OB (Figure 7) in 


Ag B 
FIGURE 7 


the direction indicated by the arrow, the muscles, wherever their at- 
tachments may be, will tend to impart a movement to the “trunk” in 
the direction indicated by the arrow. In view of the rigidity of the 
triangle OAB, the trunk will rotate around O. Only if OA is not 
rigid, and can be bent, can it push the body forward by extension. 

In case of a quadruped this difficulty does not exist. The impart- 
ing of a movement to the front extremity, as indicated by the arrow 
(Figure 6) tends to rotate the trunk counterclockwise around B , but 
this is opposed by the supporting action of the rear extremities. The 
counterclockwise motion around A , which tends to be produced by 
rotating the rear extremity as indicated by the arrow, is counteracted 
by the moment with respect to A of the force of gravity. 
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Running on extremities without middle joints. Using the nota- 
tions of Figure 6, let the extremities which are extended forward be 
uniformly accelerated through the action of forces of constant mo- 
ments in the direction indicated by the arrows. Put 


d¢/dt=$=o. (1) 


The linear velocity v of the points A and B, and therefore of any 
point of the trunk is given by 


1G) < (2) 
The acceleration is equal to 
ptt lay, (3) 


and if M is the mass of the whole system, we may approximately con- 
sider that per extremity the centrifugal force is equal to 


1 Mlw?. (4) 


Assuming that the other pair of extremities (Figure 6, broken lines) 
do not support the body during this phase of the locomotion, we find 
that the component of the force of gravity along I is 


4Mgsin¢, (5) 


where g is the acceleration constant. 
When (4) becomes greater than (5), the system detaches itself 
_ from the ground. Hence at the moment when this happens, we have: 


lo? =gsing. (6) 
Since d?¢/dt? = const = a, we have 
b = go + wot + § at?, (7) 
where ¢, and «, are the initial values of ¢ and w,fort—0. Hence 


oO — 0, al. (8) 
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Introducing (7) and (8) into (6) we find: 
1 (aw. + at)? = g sin (do + wot + $at?). (9) 


The left side of (9) is-represented by the full line in Figure 8, the 
right side by the broken line when a = 0; ¢o < 2/2, as usually will 
be the case. Then the lines intersect as shown in the Figure. In that 
case at the moment ft, the system separates from the ground. 

When the system has separated from the ground its center of 
gravity describes a parabola until an extremity again touches the 
ground. We thus have a sort of running, for running is defined as a 
locomotion, in which the animal is at times completely detached from 
the ground (Steinhausen, 1930, p. 220). Because of (3), equation 
(6) may be written: 

v?/l=gsing. (10) 

The horizontal component of the velocity of the body is always 

less than v. Hence, for the average velocity v of running 


V<v. (11) 
But from (10) it follows, that 
{lg <A e-Oreeya< LO. (12) 
Hence a fortiori 
w<clgsv< Vig. (13) 


With 1 ~ 100 cm. we find 
v < 300 cm sec". (14) 


With such a mechanism of running, a horse could not run faster than 


3 meters per second, while smaller animals would run even much 
slower. 


We may connect v with the average metabolic rate g erg. gm~ 
sec “*. While the system is still on the ground, we have: 
v=vsing;singd=v/v. (15) 
Equations (10) and (15) give 


v= v3/gl. (16) 


Neglecting the mass of the extremities as compared with that of 


the trunk, the kinetic energy which must be supplied to the system 
is approximately 


4Mv?, 
When the system detaches itself from the ground, its center of 
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gravity describes a parabola. The horizontal distance x which it cov- 


ers until it reaches the same level above ground as it had at the mo- 
ment of detachment is 


x= (2v?/g) sin dcos¢. (17) 

The time it takes to cover that distance x is equal to t = x/v. Hence 
from (17) and (15): 

t=2v cos ¢/g. (18) 


We may consider ¢ as the duration of each “running step”. If all the 
kinetic energy 4 M v? is lost after each step, due to impact with the 
ground (which is the most unfavorable assumption), then the amount 
of energy supplied per second must be 


My gMv 
Qt 4 cos o 
This must be less than Mq. Hence, cos ¢ being of the order of 1: 
gv/4=q or v=4q/g. (20) 
Relations (16) and (20) now give 


(19) 


v = 64 g*/g'l. (21) 
With g ~ 10° erg gm" sec", and 1 ~ 10? cm., we obtain v = 10? cm 
sec"!. Hence from the point of view of energy supply this mechanism 
of running is adequate to give observed velocities. It is, however, in- 
adequate from purely kinematic reasons, as we have seen above. 

We may argue that the system detaches itself from the ground 
when the vertical component of the centrifugal acceleration (v?/I) 
sin ¢ exceeds the earth’s acceleration g. This gives 


v/lsingdg=9. 
Substituting for v its expression (15), we find 
v = vVglsin ¢ 


which gives the same results as (13). 

Thus far we have neglected the masses and the moments of in- 
ertia of the extremities. In the following we shall take them into 
consideration, and therefore shall derive here some simple relations. 

Considering an extremity of length / and thickness b as a very 
oblong: ellipsoid of revolution, we find for its moment of inertia I 
around an axis passing through the center of gravity 
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4 
1 =—ap1b?(4l? + $b?) -- (22) 
120 
Let 
bV2/l=e; 2°=—el. (23) 
Then 
Cade 
I=—ap&(1+)b. (24) 
960 
Hence, putting 
4. 
a@ = — _ np &*(1 +e?) (25) 
960 
we have 
I=al’. (26) 


For similarly shaped extremities, their moments of inertia vary as 
the fifth power of their length. 
With « ~ 107, p ~ 1, we have 


aw 10¢* gmem?. (27) 


Walking on extremities without middle joints. Consider now a 
locomotion on jointless extremities in which the system does not de- 
tach itself from the ground, that is when a is so small, that t, in Fig- 
ure 8 is to the right of the maximum of the broken line. This type of 
locomotion is usually called walking. We still assume the same 
symmetry condition for the movements of the rear and front extrem- 
ities as before. During each step the trunk of mass M is lifted by the 
amount /(1 — cos a), requiring an amount of work 


Mgl(1 — cosa). (28) 


Let the extremities swing with a constant average angular velocity 
«@. For the linear velocity of the trunk we have still the equation (2). 
The kinetic energy acquired by the trunk during each step is 


4Mv? = 4MPo?. (29) 
The kinetic energy acquired by each extremity is 
I@? = al'w?. (30) 
Hence, for all four extremities it is 
Aal'q?. (31) 


Let all of that energy be lost at the end of each step. If the forward 
moving extremity suddenly becomes fastened (through friction) to 
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the ground upon contact, and the connection of the “propelling” ex- 
tremity, which at the end of the step is directed backward, with the 
ground remains unsevered, such a total loss will happen, because we 
have at the end of a step a rigid triangular structure, fastened to the 
ground. 

The total energy expenditure per step is 


Mgl(1 — cosa) + 4MPw? + 4alio?, (32) 
The distance travelled during each step is 
$— 21 sin a. (33) 


While this distance is traveled, the extremity swings through an an- 
gle 2a with a uniform angular velocity w. Hence the duration of 
each step is 

f= 20/m:. (34) 


The average velocity during the step, which is the same as the aver- 
age velocity of locomotion, is 


_v=s/t,= (lwsina)/a, (35) 


whence 
w=av/lsina. (36) 


Substituting (36) into (32), the latter becomes 


a? y2 4aa? B y2 


Mogl(1 — cosa) + 4M — : (37) 
sin? o sin? a 
The total energy supply per step is 
2Mqi sina 
Mat, = a Aa (38) 
v 
because of (34) and (36). 
Expression (37) must be less than or equal to (38). Hence 
1 Mo? v? = 4a0? | v? : ee a (39) 


Se ide 2 sin? a sin? a v 

The left side of (39), plotted as a function of |, is a cubic which be- 
comes steeper with increasing v and also is displaced upwards as 
v. increases (Figure 9,L). The right side is a straight line, the slope’ 
of which decreases with increasing v.- For sufficiently large values ° 
of’v, the two lines will not intersect, and relation (39) is impossible. 
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LM 292/2sina 


FIGURE 9 


The upper limit of v, for which the equality of (39) just holds, can 
be found in the following way. When v has its highest possible value, 
Vm, the lines L and R are tangent. Hence 
12 ao? I? v,2 2Mqsina 
Molo. 2) ee 
sin? a Um 
In order that this equation should have a real root we must have 
(2 Mq sina) /Vm — Mg(1 — cosa) > 0 

or 

Um < (2qsina)/g(1 — cosa). (40) 
For a = 2/2, (sin a)/(1 — cosa) =1. Fora — 0 it tends to o as 


2/a. Hence we can have as high a Vv», aS we wish by taking sufficient- 
_ ly small values for a. For small a’s (40) gives 


Um =4q/9 0;40=49/9 Un. (41) 


With q ~ 10° erg gm™ sec", in order to obtain v, © 10? cm sec we 
must have 


a=107 6°. (42) 
Equation (36) gives for small values of a 
owv/l, (43) 


With J ~ 10? cm we find w & 10 sec”, and therefore from (34) 
t, © 107 sec. 


The frequency of steps should be about 100 per second, while the 
length of the step is found from (33) to be s = 2la ~ 10 cm. With 
L~ 10 cm we have w & 10? sec”; f, © 10° sec, s © 1 cm. 

Contrary to the process of running, discussed above, walking on 
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jointless extremities is not limited by any kinematic considerations. 
It is not limited by considerations of energy supply, either. However, 
in order to obtain plausible values of the velocities of locomotion, it 
is necessary to have such high frequencies of steps, that two new lim- 
iting factors will enter. One is the inertia of the extremity. With 
frequencies of the order of 100 sec”, the acceleration will be so high, 
that normal muscular forces will not be able to produce them, unless 
we limit ourselves to very small organisms with small masses, like 
insects, to which the above considerations do not apply anyway. An- 
other limiting factor is the muscular viscosity (A. V. Hill, 1927) ; 
which again will prevent displacements of large muscular masses 
with such frequencies. 

All in all, the above considerations indicate the reasons of in- 
adequacy of jointless extremities for quadruped locomotion, and 
throw some light on the reasons why animals with such extremities 
do not occur. 


Running on extremities with one or more middle joints. Now 
let us consider a different type of running, when the extremities 
have one or more middle joints. By folding and unfolding the ex- 
tremities, the animal may repel the body from the ground in any def- 
inite direction. The necessity of considering four extremities, dis- 
cussed on page 18, now does not exist, and we may simplify our con- 
siderations, by considering the dynamics of two extremities only. The 
locomotion of a quadruped may now be considered as a locomotion 
of “two bipeds connected in series’”’. 

Let the repelling action of an extremity impart to the center of 
gravity of the body a velocity v, forming an angle ¢ with the hori- 
zontal, and cause the system to be detached from the ground. The 
center of gravity then describes a parabola, and the horizontal dis- 
tance traveled by the center of gravity until it reaches the same hori- 
zontal level is 


X, = (2 v?/g) sindcos¢. (44) 
The time it takes to cover that distance is 
tp = %,/v cos d = (2 vsin $)/g. (45) 


If we neglect for the time being, the length of the extremity, then f, 
is essentially the duration of the step, while 


V=VCOS ¢, (46) 


is the speed of running. 
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The energy available per step is 
Mq t. = (2 Mqv tan $)/g9. (47) 


The total energy imparted to the body, neglecting the work of the 
necessary lift during the “unfolding” of the extremity, is 4Mv?. Not 
all of the energy is lost at the end of the step. If the forward swing- 
ing extremity hits the ground at an angle 2/2, then approximately 
only the vertical component of the velocity v is lost. The amount of 
energy lost is then, because of (46): 


1M v? sin? 6 = 4M v? tan? ¢. (48) 
This must be equal to or less than (47). Hence 
v S4Aq/gtan¢. (49) 


If with g ~ 10° erg gm™ sec we wish to have v ~ 10° cm sec, we 
must take tan ¢ © 0-3 or ¢ © 15°. It is noteworthy, that examining 
photographic pictures of galloping horses, cats, dogs, etc. (Muy- 
bridge, 1887), we find in many cases, that during the “flight phase’, 
when the animal is completely detached from the ground, the angle 
formed by the extremities with the ground is of that order of magni- 
tude. 

Several improvements can be made in the above reasoning. First 
of all we should bring into consideration the length of the extremities, 
which does not enter into equation (49). We shall subsequently in- 
troduce the length of the extremities by considering the muscular 
force necessary to impart motion to the body. But quite apart from 
any such considerations, the length of the extremities may be intro- 
duced in the following manner. . 

Let the forward swinging extremity hit the ground in such a 
position, that the line OA (Figure 10) is vertical, so that approxi- 


FIGURE 10 


mately only the vertical component v sin ¢ of the velocity is lost. 


Then the horizontal velocity ¥ remains the same, and as the extremity 
now unfolds, it communicates to the body the lost vertical momen- 
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tum, which is regained when the extremity is in the “unfolded” posi- 
tion AO’. During this period the body moved the distance OO’ 
l cos ¢, and this takes a time 


t, =Ilcos $/v cospg=l/v. (50) 
Hence the total duration of the “step” is, because of (45): 
2vusin¢ 
t). + ¢, = ———__+ -, (51) 
g v 


and the energy available per step is now 


ea ee Se ee an Mal cos¢ 


(52) 
g Vv g Vv 
Setting this equal to or greater than (48), we now have: 
= 2qvtan Leos 
kv ee ae ee (53) 


g Vv 


which, for the case of equality is a cubic equation in v, and expresses v 
as a function of 1. The value of v obtained from (53) is higher than 
the value obtained from (49). That value increases with increasing | . 

We have neglected so far the mass of the extremities as com- 
pared with that of the body. In most cases, except in man, this is 
probably justified for a rough approximation. 

If the extremity which swings forward hits the ground not when 
perpendicular to it, then part of the horizontal component v cos ¢ of 
the velocity will also be lost. If we assume that all of the kinetic 
energy is lost at the end of each step, then we have, by requiring now 
that (47) be greater than or equal to 4M”: 


$Mv? S (2 Mqv tan ¢)/9, (54) 
or 


4 q tan ¢ cos? ¢ 
g 


IIA 


v - (55) 
which gives the same order of magnitude as (49). — 

The work of lifting the body during the process of unfolding the 
extremity is also in general of the same order of magnitude as the 
kinetic energy 4Mv? and its introduction does not alter the situation 
‘appreciably. 

.-. he fundamental difference between 1 this mechanism of running 
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and the one discussed above for rigid extremities, is that in the for- 
mer case the duration of the step decreases with increasing v , while 


in the latter case that duration increases with v. This is seen for the 
former case by combining (34) and (36), which gives 


t,= (2lcos¢$)/v. 


For the latter it follows from (45). Thus with folding extremities, 
the higher v , the higher the amount of energy available per step. 


Muscular forces in running. All the above discussions apply to 
continuous running and take into consideration only the continuously 
available supply of energy. Now we shall consider the problem from 
another angle, namely from the point of view of the mechanical forces 
available to produce the necessary accelerations. 

First consider the case of walking on rigid extremities, whose 
mass we again shall neglect. The length of each step is given by (33), 
while the duration is equal to t; = (21 sin a) /v. The angular velocity 
is given by (36). Hence the average total linear velocity is 


v=lw=a?/sina. 

Since that average velocity is reached approximately within a 

time 4 ¢, , beginning with zero velocity, the average acceleration is 
v= 2v/t) =a v?/l sin? o. (56) 

Hence the average force is 

M «v?/l sin? a. (57) 
This force is supplied by the muscles of the extremity and is approxi- 
mately proportional to their cross-section, other conditions being 
equal (A. V. Hill, 1927). The muscular force is opposed by the weight 


of the body which is equal to Mg. Denoting by s the average thick- 


ness of the extremity and by f a coefficient of proportionality, the to- 
tal force is equal to . 


fs? — Mg. _ (58) 
Equating (57) to (58) we obtain 


ee (fs? — Mg)l ms Ly 
pe sin a. (59) 


The force exerted by a muscle is of the order of 10 Kg em & 10° 
dyne cm. Hence f ~ 10’ dyn cm. With sinaw 1 ,ar~1;3sw 10 
cm, M ~ 10° gm; 1 ~ 10? cm; we find Mg ~ 0-1 fs?. Thus the Mg 
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term may be neglected, and equation (59) gives v ~ 10° cm sec. 

A sustained velocity of that magnitude will not be possible, how- 
ever, because for g ~ 10° erg gm- sec we must have much smaller 
values of a and sin a, as we have seen in discussing expression (41). 

In the case of running with non-rigid extremities we have for 
the duration of the action of the force the value t, given by (50). 
Since, to a first approximation, only the vertical component of v is 
affected in that case, we see that during the time 1/v the veloc- 
ity changes from 0 to v sin ¢. Hence, the average acceleration is 
v? sin ¢/l and the average v force is Mv? sin ¢/l. Putting this equal to 
fs? — Mg, and remembering (46) we find 


=: (fs? — Mg)l 
1= 


OS ¢, 60 
M sin ¢ icy 


which is similar to (59). For ¢ © 15°, which, as we have seen, (p. 

26) is needed for a sustained velocity, we find cos ¢/\/sin ¢ © 17, 

and thus with non-rigid extremities velocities of 10? cm sect may be 

sustained continuously both from the point of view of energy supply 

as well as from the point of view of muscular forces available. 
Equation (60) may be written 


Gc Sass ; 
= COS ¢. (61) 
M sin oe sin ¢ 


The quantity s*/ is the volume of the extremity, which is roughly pro- 
portional to its mass. Hence for similarly built animals s7l/m = y = 
const. If the second term is small, as compared with the first, which 
as we have seen is the case, then 


sa /o COS ¢. (62) 
sin ¢ 
Hence, for similarly built animals the maximum velocity is inde- 
pendent of the size of the animal, other physicochemical conditions 
(which enter in f) being the same. This result has been previously 
obtained in a different way by A. V. Hill (1927). 

On the other hand, if we assume, as in section JII, that s* is pro- 
portional to M , so that, denoting by k a constant s? = kM, then, again 
neglecting the term in g: 


fki 


; (63) 
sin 


_—— 
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The velocity now is independent on the mass. of the animal, but 
proportional to the square root of the length of the extremities. Simi- 
lar relations have been derived previously from purely dimensional 
considerations involving some doubtful assumptions (Thompson, 
1942). A comparison of such a relation with observations presents 
great difficulties in view of the scarcity of data. Hitherto almost the 
only available data on the speed of animals are given by the series of 
instantaneous photographic pictures published by E. Muybridge 
(1887). They give interesting insight into different modes of loco- 
motion, such as walking, trotting, galloping. In no case do we have 
the certainty that the maximum speed has been observed in any of 
the locomotions. For comparison with equation (63) we must choose 
the speediest form of locomotion—galloping. A few data could thus 
be found in Muybridge’s publication, and they are represented on 
Figure 11. For a greyhound the observed velocity is however much 
too high. 
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FIGuRE 11. Velocity v of galloping plotted against length J of extremities on 
logarithmic scale. The straight line has the theoretically required slope 1/2. 
ene 2 — goat; 3 — elk; 4 — horse; 5 — camel. Data from Muybridge 


It would certainly be of interest to obtain more detailed data on 
various animals by means of a cinema film. The quantity f could in 
principle be measured on each animal and the cross-section s? can also 
be directly determined. Thus it would be possible to directly compare 
equation (60) with observations. 

Before any worth-while discussion of the available data can be 
made, the theory will have to be considerably improved. Instead of 
using average velocities and accelerations, it will be necessary to es- 
tablish more exact, although still approximate, equations of motion 
for the folding and unfolding of an extremity. The relative dimen- 
sions of each part of the extremity, in other words, the location of 
the middle joint, may then be determined by the requirement of a 
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maximum possible velocity, considered as a function of those dimen- 
Sions. The next extension would be to determine from similar con- 
siderations the optimum number of joints. Besides considerations of 
locomotion, other considerations will enter into that problem. An ex- 
tremity ending in a thin stump would be of little use on a soft ground. 
This alone, regardless of other reasons may require an extra joint, 
forming the foot. 

A more detailed study of the losses of the energy due to impact 
against the ground, will also have to be studied. This also will con- 
tribute to the determination of the structure, and hence of the shape 
of the extremities. An approach in that direction is made in section 
V. 

The requirements of maximum short time velocity or of maxi- 
mum acceleration are different from those of sustained running. The 
study of this aspect of locomotion will lead into a theory of jumping 
and require considerations of extremities of unequal length. Even 
the role of an additional extremity in jumping—the tail, may be 
brought quantitatively into those studies. All such studies will deter- 
mine more closely the external shape of the animal. But even with 
the present crude approximation we have the following result. 

For a given mass of the animal, the size and shape of its trunk 
are determined by relations (7) of section III, while the thickness of 
the extremities is determined by relation (8) of that section, the 
length of the extremities is given by equation (60) of this section. 
The general outside shape of the animal is thus determined quantita- 
tively to quite a large extent even by those simple considerations. 

In addition to that it may readily be seen that the total length 
of the neck must be of the order of magnitude of the length of the 
extremities, in order that the animal could reach food on the ground. 
The thickness of the neck will be determined for a given length, by 
the weight of the head. The latter in its turn is approximately pro- 
portional to the weight of the brain. It is empirically known (Crile, 
1941, p. 210), that the brain weight of animals is approximately pro- 
portional to the total metabolism*) gM. Hence, q and M determine, 
together with the shape of the trunk and the thickness and length of 
extremities also, the size of the head and the length and width of the 
neck. 

With the exception of the tail, we thus have a fairly complete 
description of the general shape of a quadruped in quantitative terms. 

The approximate proportionality of the weight of the brain and 
of gM still awaits a theoretical explanation. Perhaps it may be found 


* Crile does not give the source of his data for q, and we feel therefore, that 
the above relation must be taken with due reservation. 
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through the consideration that the more complex a brain, the larger 
the number of innervated muscle fibers and other end organ units, the 
tonic activity of which consumes energy. 

The above mentioned linear relation between gM and the weight 
of the brain makes it rather pointless to assume, as we have done pre- 
viously. (Rashevsky, 1943a), that the brain weight is to be considered 
as an independent determining parameter. It also makes untenable 
the previously suggested hypothesis, that the degree of nervous co- 
ordination available is responsible for the use of two rather than four 
extremities for walking and running (Rashevsky, 1943a). The rea- 
sons why most animals use four extremities for locomotion, when 
two seem to be sufficient, must be sought elsewhere, and may well lie 
in some details, as yet not studied, of the dynamics of locomotion. 

As regards arthropods the larger number of extremities in those 
animals will likely be explained by the “relative roughness” of the 
ground (Rashevsky, 1943a; p. 40) and other similar considerations. 
The locomotion of jumping insects, like fleas or grasshoppers is more 
aptly described as biped. 

All the above holds on the assumption of the constancy of other 
parameters, in particular of f. The latter is largely determined by 
the physicochemical constituents of the animal. It may well happen 
that f may itself be considered as an independent determining para- 
meter of the organism. This does not alter anything in the general 
situation. We determine the shape of the organism in terms of con- 
tinuously varying biological parameters, which describe the physico- 
chemical constitution and the physiological functions of the organ- 
isms. 


V. Loss of Energy Due to Impact of Extremity Against the Ground 


In this section we shall study a few of the simplest cases of the 
problem. Further extensions of this investigation will be evident. 

Let us compute the loss of energy, due to impact against the 
ground of an extremity, without knee-joint, directed forward. Let 
us first neglect the mass of the extremity. Let (Figure 12) the center 
of the gravity of the whole system be at S, , the middle of the line 


be} (b, 0 


aes 


A 


Ficure 12 
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OO’ , of length 2L.. Let the coordinates of S, be x and y , and the ve- 
locity components just before the impact be 


After the impact, which is supposed to take place simultaneously 
with the front and rear extremity, so that OA || O'A' and OO’ is hori- 
zontal, we have, denoting by 2, the abscissa of the point of impact A’, 
which after the impact remains fixed: 


(eel ee egy fe = (0), (2) 
or, putting 
OP et B aaa 
(x —2)?+y2?-—P=0. (3) 


We shall apply the Lagrangian equation for the calculation of 
values of x and y after the impact (Appell, 1924, p. 531). To this 
end we introduce new variables 

i=“; q%=(&—x)?+y¥—-P. (4) 


Due to the constraint (3) introduced by the impact, we have at the 
moment of impact and at all times thereafter 


Q2— 0. (5) 
Denote by T the kinetic energy of the system 
T=4M (a? + y’). (6) 


By making the substitution (4), T becomes a function of qi, @, on 


and qg.. We now have, denoting by the superscript 0 the values 
just before the impact, and by the asterisk the values just after the 


impact: 
or \* oT \° | 
egies oes) = 0) (7) 
ads 0 

From (4) we have 


y=qi+P— (z—x)*. (8) 

Differentiating, dividing by 2y, squaring and again using (8), we 
find: 

“S Q2"Qz? Ec (x 197)" qn? ae 24242 (Zo - 1) 

> se ac = 


= (9) 
q2” +?— (x ae Gaye 
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Introducing (9) into (6), we find, remembering that, because of 
(4), fe ad h ? 


| ofp Cee ae ee _ (10) 
go? + Ss Tes qi)? 
Hence 
oT . M(x my Qs ae M242 (ao — q;) 
aq, Ga? (eda)? 


To obtain (oT /eq:)* and (@T/2q:)° we must substitute into (11) the 
values which q,, 91, 2 and qd. have just after and before the impact, 
correspondingly. We shall denote those values with an asterisk and 
a superscript ° respectively. But q.*? is zero, because go* becomes zero 
at the moment of impact and remains zero thereafter. Keeping this 
in mind, as well as the fact that the denominator of (11) does not 
change during the impact and is equal to /? sin? ¢., where ¢» is the 
value of ¢ at the moment of impact, we have: 


olen OT Ne c i 
Shi eae | oa M (q.* — 1°) 
0g: "1 


Ma ia) (q* ae 1°) oe Mq.°Qq. (x oy a!) 
I? sin? do 


(12) 


Now we may substitute into (12) the original variables x and y. We 
notice that 


a." — 1° =a" — 2° =v" —v°; (4% — a1) =L eos go. (13) 


‘ The value of qo°qs° just before the impact is otained from (4). 
Differentiating with respect to t, we find that before the impact: 


G24 = yy — (@— 2) a. (14) 
but at the moment just before the impact,  — x = 1 cos &, y = 
lsingd.,x=v,y=—w. Hence 

G2° G2" =— 1(v cos bo + wsin do). ~ (15) 


Equation (15) shows, that, since at the moment of impact Qo° 


becomes zero, g2° becomes negatively infinite. This is: obvious, since 
the lel, changes discontinuously. 
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Introducing (15) and (13) into (12), and the latter into (7) we 
find, after rearranging: 


v* — Vv = — COS go(V COS do + win do) <0. (16) 


The expression in parentheses is the projection b; on J of the to- 
tal velocity » of the center of gravity. Hence the loss of the horizontal 
component of b is equal to the horizontal component of »; . 

Solving (3) for y? and differentiating with respect to ¢, we find, 
remembering (13): 


y= ‘a *= z* ctn do. (17) 


To compute the total loss in kinetic energy due to the impact, we 
remember that 


T*=1M (x + y”); T°=1M (v2? + w). (18) 
Hence 


T° —T*=4M (v? + w? — 2? — yy”). (19) 


Solving equation (16) for v*, introducing it into (17) and then 
introducing v* and (17) into (19), we obtain after rearrangements, 


AT = T° — T*=4M (v cos do + WSIN do)? = $M,” . (20) 


Thus the total loss of kinetic energy is equal to the kinetic energy 
of a system, having a velocity equal to the projection on | of the actual 
velocity of the system. _ 

AT becomes equal to 4M (v? + w?) = T°, that is all kinetic ener- 
gy is lost, if » coincides with t; , that is if 


Ww 
= = tan go5 w= 0 tan go. (21) 
AT is zero, when 
oy) 
tan Por a (22) 
w 


that is, when » is perpendicular to /. 
If the center of gravity moves before the impact only under the 
influence of gravity, having at the moment t = 0 the coordinates 


x2 = 2, and y = h and the velocity component z= vand y = 0, thus 
describing a parabola, then we have 


Z=vt; y=h— tot. (23). 


At the moment of impact t , y=J sin ¢.. Hence 
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h — 4gt.? =1 sin do 


or 
2h — 2l sin 
ty = J aoe Pees (24) 
g 
The velocity y is equal to — gt, and therefore at ¢, it is equal to 
— w= — gt =— 29 (h — Isin go) (< 0). (25) 


Introducting this into (20) we find 


AT =14M (v cos 6. + V2g (h — Usin $0) SiN go) « (26) 


In particular we may consider the following case. Let the sys- 
tem be propelled to the right by the rotation of the extremities clock- 
wise. Consider that at the moment when ¢ = 2/2, and S ‘is at its 
highest position, the lower end of the extremity becomes disconnected 
with the ground, sliding on the latter without friction. Then S will 
move on a parabola, its coordinates being given by (23). The ex- 
tremity, paired to the propelling one, is swinging forward. We as- 
sume that it also drags on the ground without friction, until its angle 
¢ becomes equal to ¢).. Then it suddenly becomes attached to the 
ground. A similar situation will arise, if, in order to avoid rubbing 
the ground with the forward swinging extremity, the body is slightly 
swayed to the side opposite to that extremity. When ¢ = ¢,, then the 
body weight is applied to it in full. 

We now have h = 1, and (26) gives 


AT = 3M (v cos ¢o + V2i g(1 — sin go) sin go). (27) 


AT = 0 for ¢ =2/2. The loss AT is the less, the smaller the “steps”. 
If we take into consideration the mass of the extremities, the ex- 
pressions become more complicated. 


VI. Suggestions for an Approximation Method for Solving the 
Equations of Motion of a Chain of Linked Levers 


The equations in question are of the form (O. Fischer, 1906; 
Winkelman and Grammel, 1927) 
obi +S Oy (be + ox?) f (oi — ox) =U(¢i) , (1) 
k 


where ¢; are the angles determining the positions of the individual 


levers of the chain, f and u are continuous functions, and a, are con- 
stants. 


Let during the interval 0 — t, the average values of $i be i. 
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Then the average values of ¢; at the moment t, are : = $4; ty. The 


average accelerations are ¢ = At Substituting those averages into 


(1), we have 
Oy b/ty +S (Gu /ty + bx? )f Adits — Abets) =Ul(dit,). (2) 
k 


These are algebraic equations, which determine the values bi as 
functions of t,. 


If at t= 0, 6: = $i 2 0, then we put fort=—t,: 


re boi + bi - bi — bos 
pen ¢ oe 2 bes. (3) 


2 be 


Substituting (3) into (1) we again obtain algebraic equations, ex- 


pressing ¢; in terms of t, and $v: - 
This method may work if the velocities and accelerations are 
monotonic functions of ¢ in the interval 0, ¢, and if the accelerations 


do not vary too rapidly in that interval. 


VII. Flight of Birds and Insects in Relation to Their Form 

Even a remotely exact theory of flight of birds and insects meets 
with almost insuperable difficulties at present. In fact the relatively 
simpler problems of aircraft design are still based on semi-empirical 
foundations, since a complete synthetic theory of aerodynamics is still 
wanting. The relative simplification of the problems of aircraft en- 
gineering is due partly to the circumstance, that the two principal 
types of aerodynamical structures involved are the ordinary wing 
and the propeller. The action of both involves to a large extent sta- 
tionary aerodynamical states, when continuous operation is consid- 
ered. The situation with birds is different. While the wing repre- 
sents an aerodynamical profile, yet its movement through air is any- 
thing but stationary. Inertial forces due to acceleration, which do 
not enter into the theory of the airplane wing, probably play here a 
very important role. 

From studies by Marey (1893) and from pictures by E. Muy- 
bridge (1887) one gets the following general impression of the me- 


chanism of the flight of a bird. 
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During the downstroke, (Figure 13a) the wing moves down and 
forward, with the velocity v. The lift A has thus a component di- 
rected forward (A,). The resistance W has a horizontal component 
W, directed backward. Apparently A, — W.>0, and causes the bird 
to move forward. The vertical components A, + W, compensate at 


FIGURES 13a AND 136 


least partly for the force of gravity Mg , M being the mass of the bird, 
g the acceleration of gravity. Actually A, + W, is less than Mg , be- 
cause, according to Marey, during the downstroke the bird looses 
slightly in altitude. 

During the upstroke the wings move backwards and upwards 
much more slowly than during the downstroke. The bird moves for- 
ward, due to the momentum imparted previously by A; — W,, with 
a velocity v, (Figure 13b), which gradually decreases due to resist- 
ance. Because of a larger angle of attack a, the lift A, , which is now 
directed vertically, is greater than A, before, so that A, > Mg > A,, 
and the bird again gains altitude. 

While in a sense the wing acts as an aerodynamic profile both 
during the downstroke and during the upward stroke, yet in no case 
can the ordinary theory of an airplane wing be applied with any ac- 
curacy. For neither the velocity v in Figure 13a, nor v, in Figure 13b 
are uniform. Nor is the direction of the velocity with respect to the 
wing constant. In fact in both cases the tip of the wing has a much 
higher velocity than the point of attachment which remains station- 
ary with respect to the body of the bird. This suggests that the wing 
may be acting as a propeller blade (Pettigrew, 1875; cited from M. 
H. Fischer, 1930, p. 340), rotating approximately around the axis’ 
OO’ (Figure 18a) and propelling itself in that direction. The verti- 
cal component then provides for the maintenance of altitude, the 
horizontal—for the forward movement. But again the fundamental 
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difference between a theory of such an action and that of a propeller 
is that the former does not correspond to a stationary state of rota- 
tion. Moreover the two wings must be considered as two propeller 
blades moving in opposite directions and their mutual disturbance will 
be very complex. Thus we have in a bird a new type of aerodynamic 
lifting surface, which is neither a stationary moving wing nor a pro- 
peller, but a system sui generis, the theory of which must still be 
developed. 

From Marey’s studies it appears that the flight of insects is 
based on a somewhat similar mechanism. 

At the end we shall suggest a rough, but somewhat detailed 
analysis of the above discussed process of flight. Now we shall con- 
fine ourselves to very rough and general approximations, based 
largely on dimensional considerations. 

Whatever the details of the actual mechanism may be, we may 
consider that the total velocity of locomotion v of the animal through 
air is less than, but of the same order of magnitude as the velocity 
v of the wing through air. Thus 


VD, (1) 


Denoting by L the length of a wing and by S its width, by p the 
density of air (p = 1-2 X 10°? gm cm‘) and by ¢, and «@, the lift and 
resistance coefficients of the wing, we find for the force of resistance, 
F,, approximately 

F,=4 Cy pLS® . 1 (2) 
If the frequency of wing movements is », then the average velocity . 
v is 
yvnwlLy. (3) 
Since this velocity v is acquired approximately during a time 1/» , the 
average acceleration is 


voLy, oe (4) 


If we put the density of the wing equal to p’, and its average thick- 
ness to A, we have for its mass M,,: 


My © LSp'A. (5) 
Hence the average inertial force is 
F,=M,0% L?S »? Ap’. (6) 


The total force F; + F, must be exercised by the body in order 
to move the wing with the desired velocity. This muscular force may 
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be put as proportional to M*/*, since it is proportional to the cross 
section of the muscles. Hence we have from (2) and (6), denoting 
by a’, and f two coefficients: 


W's L?SAy’p’ is 4 CypLS We @ f yy ae (7) 
or, because of (3), with a, as a new coefficient and remembering that 
there are two wings, so that F, must be doubled: 

Ay SA p. (Va Cup LS VEE f M?% : (8) 


The lifting force of the two wings is equal to ¢,pLS v* and this 
must be equal to Mg. Hence 


Mg 
LS « : (9) 
Cu p UV" 
Solving (8) and (9) for L and S, we find: with e = ¢,/¢,: 
a,MAgp AA 
Li« ee (10) 
Cap (f{M ane g) het eg 
wy M 
and 
fM?" — «Mg 
Ss e————_ (11) 
ay Ap’ vy 


A relation, determined by questions of mechanical strength, ex- 
ists between L,S, A and M. The thickness A must be the greater, 
the greater L and M. The situation is similar to the problem of the 
thickness of a branch of a tree or the length of the trunk of an ani- 
mal. Thus 

E(G5S5 435A). = 08 (12) 


Equations (10), (11) and (12) give us the size of the wing in terms 
of M and v., since v « v. 

Equation (10) shows that there is a limit of the mass of a bird. 
For when M > f?/e? g?, L becomes infinite. The actual limit is of 
course much below (f/eg)*, since for structural reasons the length 
of the wing cannot exceed a certain size, which is of the order of mag- 
nitude of the size of the body. With f ~ 10°, which we obtain by con- 
sidering that the cross section of the muscles in question is about 1% 
of the cross section of the body, (f/eg)* © 10° gm = 10° kg. 

If d « L, then L and S cannot be determined from equations 
(8) and (9), since only the product LS-enters into them. In that 
case equation (9) gives us LS, while equation (10) establishes a re- 
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lation between the different parameters, which must hold, in order 
that flying should be possible at all. In particular, putting 4 = pL, 
we have — : 2 


aig B p' 


p eat 
wm” 


which shows that the larger M the greater must be c,. Hence for 
larger animals, like birds, the profile of the wing must be much more 
perfect an aerodynamic profile, than for insects. Actually this is the 
case, the wing of a bird being much more perfect aerodynamically, 
than the wing of an insect. 

Actually in equation (8) we may have more generally a “less 
than” sign instead of an equality sign. This, with 4 = BL will result 
in a “greater than” sign in expression (13). Thus the conditions 
imply inequalities rather than rigid equalities. 

To give one illustration of the above: The load per square cm in 
birds is of the order of 0-4 — 1-5 gm cm~ or 4 X 10? — 15 X 10? dyn 
em’. With velocities of the order of 10° cm sec’, this gives values 
of c, between 0-6 and 2. For insects, the load per cm? is “ to — 
of the above. With velocities of the same order of magnitude, this 
gives a correspondingly smaller ¢,. 

For approximately similar animals the mass varies as the cube 
of the linear dimensions |, while the area of the wings varies as [?. 
Hence the load per unit area varies as 1. The velocity v varies as ly , 
hence the resistance as /?v? = Ity?. The acceleration varies as ly’, 
and the mass of the wing as /?,, hence the inertial force varies as I*y” 
also. The sum of the resistance and inertial force must be propor- 
tional to the muscular force, which being proportional to the cross 
section of the muscles varies as /?. Hence 

PE«P, 
or remembering that l is prc ogrtional to the load m per cm? of the 
wing 


Cai (13) 


ya l/l; y « 1/m. 7 (14) 


These relations hold, however, only for similarly built animals. In 
approximately similarly built insects it is actually observed that the 
larger the animal the smaller », so that at least qualitatively rela- 
tion (14) is verified. (v. Buddenbrock, 1930, p. 353). 

If, however, we consider animals which have approximately the 
same mass but wings of different sizes we obtain a different rela- 
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tion. The inertial and resistance forces are again proportional to 
I472. But the muscular force is now roughly constant. Hence »*/* = 
Const and » «1/l?. But for constant mass, /? is inverse proportional 
to the load m per cm? of the wing. Hence instead of the second rela- 


tion in (14), we have: 
yam. (15) 


Relation (15) is very well verified by observation for some insects 
(Figure 14), although they appreciably differ from each other with 
respect to their mass. 
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Figure 14. Frequency v of wing beats of some insects plotted against the 
inverse of load per cm? of wing. Experimental points compared with the theoreti- 
cal curve. Adapted from Buddenbrock (1980, p: 354). 


It is remarked, however, that larger insects fly more rapidly than 
smaller ones. This at first seems to contradict the above considera- 
tions, since v « ly = const. But we must remember that v differs 
from 7. The coefficient of proportionality between v and v, appar- 
ently varies with size. This is still a problem to be studied. 

If we compare the empirical values of c, and c,. for aerodynam- 
ical profiles and for plane plates, we find (Betz, 1927, pp. 285-287), 
that a value of ¢, © 1 may be obtained with either. However in order 
to obtain such a high value of c, for a plate, a must be taken rather 
large. Then c, turns out to be about ten times as large for a plate as 
for a profile. For larger animals, like birds, a wing in form of a flat 
plate would require very high energy expenditure and great force, 
due to a high ¢,,. For insects c, may be about ten times lower. There- 
fore a flat plate can be used with a much smaller a , In which case cy 
will also be reduced correspondingly. Thus insects may have smaller 
¢, than birds, but approximately the same c,,. The resistance, as we 
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have seen, varies as I‘? , ¢,, determining the coefficient of proportion- 
ality. Muscular force varies as 2. The relation »y = 1/m (equation 
15) should hold over the whole range of birds and insects, if Cy» is 
about the same. 

If A « Land only LS is determined, then the determination of 
L and S separately is made in the following way. For a given LS, 
the induced resistance will be the smaller the greater L. Hence c, 
is itself a function of L and S separately. Thus, we still have two 
equations for the two variables L and S , although they will be rather 
unmanageable. In general L should be taken as large as is compat- 
ible with structural requirements, such as lead to equation (12). 

If the velocity v is to be sustained, then we have an equation 
which relates that velocity to qg. The work done by the force F,, per 
unit time is for the two wings: 


20F,=—ypLSv. (16) 
The work done by F; per unit time is approximately 
IAS 57 A pe. (17) 


Hence per unit time, because of (3): 
L?¥ SA p=Ap' Sv . (18) 


The sum of (17) and (18) must be equal to Mq, if Mq is the excess 
metabolism over the resting one. Hence 


(Ap +eypL)SvV=Mq. (19) 


Since L , S and A are determined from (10), (11) and (12), equation 
(19) gives a relation between g, M and v. Hence q and M alone de- 
termine the motion and form of the animal. 

Now we shall briefly outline a more detailed but very Fae ap- 
proach to the problem of flight of birds. 

In glancing over the graphs, representing the empirical relations 
between c, and a (angle of attack) and c, and a (M. H. Fischer, 1930, 
p. 323) we find that with some approximation, within a fairly wide 
range of a’s , we have: 


SA eb ac — CD ats (20) 


Let (Figure 15) the bird hold the wing at an angle # to the vertical 
and swing it downwards. The resulting forces will move the wing for- 
ward, so that its resulting velocity will be v , making an angle ¢ with 
the vertical. The vertical force due to the lift 4 and resistance Fy, 


is equal to 
F=F,sin¢+F,cos¢, (21) 
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FIGURE 15 


and is directed upwards. In order to move the wing down, it is nec- 
essary to apply a downward force F,. If Mw is the mass of the wing, 
and v, its instantaneous vertical velocity, then 


M,—=F,-F. (22) 
dt 
Let v, (Figure 15) be the average vertical velocity, then the du- 
ration t, of a downward stroke is approximately 
t1,=L/,. (23) 
Hence the average acceleration is 


VV /les (24) 
Averaging equation (22), we have 
Bh Ct acc 
ee , (25) 
But 
Pej fits (26) 
On the other hand for each wing: 
Fz=4o,pLS8v?; Fyox=teypLSv?. (27) 


Hence, introducing (20) into (27), the latter into (21), and the lat- 
ter with (26) into (25), we have 


M1 Vv? 
L 


= fM’?* — pLS{[A + B(6 — ¢)] sing 


+ [C+ D(6 — ¢)?] cos d}v?. = 
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Furthermore we have: 
V;=VCOS ¢. (29) 
The sum of the horizontal components of F', and F'y is: 
F’=F,cos¢—Fysing, Sen eee (30) 


and must be positive. 

If M is the mass of the bird, and R the coefficient of resistance 
of the body to the horizontal movement, we have, Se ey Vr, the 
velocity of the body in the horizontal direction 


¢ —— tt — f2 92... (31) 


Since the duration of the down stroke period is given by (23), the 
average acceleration is 


Soa L | (32) 
Equation (31) now gives 
Mv, Una ‘aa eos 
= am R Vhi e an (83) 
L nee 
We also have the requirement that, approximately at least, 
F=Mg. va (34) 


Expressing in (33) and (384) F’ and F in terms of v, LS , fh and ¢, 
we thus have four equations, (28), (29), (33) and (34), which deter- 
mine v, LS, ¢ and (<1) for a given 8. The distance traveled 
during the ee stroke period is 

= L Um 


Un ty = (35) 


V4 | 

For the upstroke we can consider the movement of the wings 

relative to the body as slow, and compute the loss of velocity and gain 
of altitude by putting 


dv 
M—=— (cy p LS + R)v%n. | (36) 

dt 
Easily integrated, (36) gives Viz as a function of time, with the initial 


conditions Vj2(0) = Vn. « Co is now to be taken for a different a. The 
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duration of the upstroke period is given by the requirement, that 
Cap LS Ving 2 Mg . (37) 


The moment when the inequality (37) ceases to hold, the new down 
stroke must begin. 

The duration t, and t. of the two periods thus known, as well as 
the average velocities during those periods, the total average velocity 


v is obtained in terms of M, L,S, a1, a etc. 
The angle a may be chosen so as to maximize v. 


VIII. The Internal Structure of Animals 

The size and shape of different internal organs may also be con- 
sidered as determined by the requirements put on them by the per- 
formance of their functions. Hence a general physicomathematical 
theory of organic form should include the internal structure of ani- 
mals as well as their external shape. A step in that direction is to be 
made in this section, which, however, must be left even sketchier than 
the preceding: ones. 

If the mass M and the maximum sustained velocity of an animal 
are given and if we denote by q* the metabolic rate necessary to sup- 
ply the energy for sustained locomotion, then, for the case of a quad- 
ruped or biped, the value of q* is given by equation (49) of section 
IV. But a continuous supply of the energy Q = Mq’* per unit time re- 
quires a definite continuous supply of oxygen. If V, denotes the am- 
plitude of the volume of the lungs, and », the frequency of respira- 
tion, then the supply of oxygen is proportional to »,V,. But the oxy- 
gen supplied goes not only for the energy supply of the running 
mechanism, but also to supply the energy Q, necessary to maintain 
the movements of the lungs, the energy Q., necessary to maintain 
the circulation of the blood, the energy Q; , necessary to maintain the 
peristaltic movements of the intestinal tract, and the energy 
QQ.) =Mq , necessary to maintain the nonmobile cells of the organism. 
Denoting therefore by a, and a., two coefficients, which are deter- 
mined by the efficiency of the lung (a,), and by the energy equivalent 
of the oxygen, we have 


a, v, Vi, =a) (Mq* + Q,+Q.+ Qi +Q). SNL 


But the maximum frequency y, is determined by V,. The latter de- 
termines the size of the musculature that contracts and expands the 
lungs. This size determines first of all the possible force acting on 
the lungs, as well as the viscous resistance of the muscles. By appro- 
priate mechanical considerations, involving possibly also the inertial 
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terms, we shall arrive at an equation 
vp —=f1(Vz). (2) 


Furthermore, if », and V, are given, the rate Q, of-energy supply 
necessary to maintain the periodic contractions and expansions of 
the lungs is also given: 

Q. =F. (11, Vz). (3) 


The adequate supply of oxygen is provided also by the circula- 
tory system. Denoting by V, the effective volume of the heart per 
beat, and by », the frequency of heart-beats, we see that »-V, is pro- 
portional to the sum of all the Q’s. Hence with a; as a constant 


veVe= a3(Mq* + Qi + Qe + Qi + Qo). (4) 


And again we have 
¥eo=Ffe(Ve), (5) 


where f, is different from f,. 

For Q. we must remember that it depends not only on », and 
V., but also on the size and structure of the blood vessels. For the 
same », V., Q- will be greater, the longer the total lengths J, of the 
circulatory circuit and the less its effective cross section s.. Thus 


Q-=F.(%» Vez bes Se) (6) 


But /. is determined by the overall length of the body L, and of the 


extremities 1. Thus 
L=f,(L,1). (7) 


Similarly, the value of s, is determined by the total size of the 
animal. For an adequate supply of all tissues with oxygen and ade- 
quate removal of waste products, there must be so many capillaries 
per unit volume, and they must be of a given size, to have the neces- 
sary surface. Hence 

Se =f2(M ya (8) 

Denoting by J; and s; the length and average radius of the intes- 
tinal tract, we find that the total area 2zl,s, must be proportional to 
(Q, + Q', + Qe + Q'i), where Q',, Q’., Qi denote the necessary 
energy supply to the lungs, heart and intestines in a resting organ- 
ism, that is when g* —0. The reason for the proportionality is the 
necessity of a given rate of resorption, which we assume to be pro- 
portional to the surface. Hence 


2al;s; = A4(Qo am QL + Q. + Q’:). (9) 


Similarly, the volume of the intestinal tract is determined by the 


48 STUDIES IN ORGANIC FORMS 


same quantity, since the latter determines the amount of food that 
must daily be held in the intestinal tract: 


Ls2@=f;(Qo+ Q14+Qe2 +03). (10) 

Also we have 
Q'.=Fi(l,, 3). (11) 
The length of the body may be determined by mechanical con- 


siderations from its mass M, as has been done in section IJI. We 
have in general 


Ler (12) 
For the thickness s. of the legs we have, as in section III: 
s?=a;M , (13) 


and finally we have expression (62) of section JV, giving / in terms 
of the maximum velocity v: 


=a = (14) 


For an animal at rest, g* = 0, and Q, = Q'.,Q-=—=Q'.,Qi=Qii.- 
The latter remains also for g* > 0, when the animal runs. 

Putting g* = 0 into equation (1) and (4) we have 14 equations 
for the determination of the 14 quantities: 


WE Vin Qt, Qc, Qi, ve; Vous bey Oat i 3 Sines Cant (15) 


in terms of M, g* and v. 

Returning now to the original equations with q* > 0 and when 
it is given by equation (49) of section IV, equations (1), (2), (8), 
(4), (5) and (6) give us 6 equations for the determination of the 6 
quantities: 


Pian Qin Gas eeu as (16) 


Thus, M and q, determine not only the general shape of the ani- 
mal, but also its internal structure. 

Strictly speaking, relations (2) and (5) will involve m, and m,, 
the masses of the lungs and of the heart correspondingly. These may 
be left as undetermined parameters. Then, according to the principle 
of maximum simplicity, we shall take for m,; and m, the smallest val- 
ues, which permit real positive solutions (15) and (16) of our equa- 
tions. That such minimum values of m, and m,, different from zero, 
must exist, is physically obvious. 

Kstablishment of the form of f, (equation 2), F, (equation 3); 
f- (equation 5), F. (equation 6), etc. form the main part of our prob- 
lem. This is to be done through a theory of lung and heart activity. 
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We may introduce equations for other organs. Thus the mass 7 
of the kidney will be determined by the maximum metabolism SQ, 


which determines the necessary excretion. The mass of the liver will 
be at least partially determined by the amount of glycogen stored, 
which may be mobilized during running, to supply the necessary q*. 


Notice that } Q', represents fundamentally the basal metabol- 
k 
ism. It may be possible that expressed as a function of M, > Q’; 
k 


will vary approximately as M?/’. , in spite of the fact that Q, « M. If 
this were so, we would have Rubner’s law derived from the above 
equations. 


In the above considerations we have entirely disregarded the 
skeleton as an internal structure. We may introduce it in the follow- 
ing manner. 


Instead of (12) we shall have, Taree by L, and s, the length 
and thickness of the spine, 


be FM se (17) 
this may be subject to the condition that the ratio 
Ls 5,337 
R= =p (18) 


would have a minimum. Expressing L, in terms of M and s, from 
(17), we find: 

R= R(M,§s8,) 
and determine s, from 


—=0. (19) 


This would be in line with the principle of maximum simplicity, 
for it will give us the minimum weight (« volume) of the skeleton. 

Relation (13) should, strictly speaking, determine the cross sec- 
tion of the bones of the extremity. The muscular cross-section may 
be argued to be proportional to the latter, because the stronger the 
bones of the legs, the stronger muscular forces and torques they may 
support without breaking. 

Equations (49) and (62) of section IV both contain the angle ¢. 
The requirement that both maximum and sustained values of v should 
be the same eliminates ¢, giving a relation between q , v and 1, which 
takes place of equation (14). 
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If we add 
M=>M,+M.+Mi+M, 


we may perhaps eliminate one of the three determining parameters. 

As a matter of illustration only, without any claim to any corre- 
spondence to real situation, we shall briefly outline a possible theory 
of the lungs and heart. 

Let us consider the plausible case, that Q,, Q- and Q; are negli- 
gible compared to Mq* + Q., and let the latter roughly follow Rub- 
ner’s law, being proportional to the surface of the animal, or, for 
similar animals, to M??. Then 


Verna Mee: (20) 


Let us also assume, that the volume V, is proportional to the mass 
M, of the lungs. The linear dimensions are then proportional to 
M,2" , and the acceleration in the periodic breathing movements are 
proportional to M,1* »,;2. Hence the inertial forces are proportional 
to M,‘ »,;2. Since the muscular forces of the lungs vary approxi- 
mately like the cross section of the muscles, they will be roughly pro- 
portional for similarly built lungs to M,?. As the lung pumps air 
back and forth, it has to overcome the resistance of the air flow 
through all the parts of the lung. That force of resistance may be 
taken as roughly proportional to the total pressure, which is propor- 
tional to the volume velocity V, », <« M, v,. The force of resistance 
will also be proportional to the total area S in which that pressure 
is exerted. If the size and number of alveoles per cm? is roughly con- 
stant, then S « V, « M,, and the resistance force is proportional 
to M,? »,. Since the muscular force must be equal to the sum of the 
inertial forces and of the forces of air resistance, we have, denoting 
by a and b two coefficients: 


M4? »,;2 + aM ,2 ve — 0M 74 — 0 (21) 


From (20) and (21) we have, remembering that V, « M Lt 
M+4/9 y,/3 ae alMs/2 ae, by,1/8 = 0 5 (22) 


The graph of equation (22) is represented on Figure 16. From it we 
see that it has two positive roots, one increasing with M , the other 
decreasing. It follows from (20) that the smaller root, which in- 
creases with M, gives a larger value of M,. By the principle of 
maximum simplicity of structure, we must therefore choose the larg- 
er root, which gives a smaller size of lung. 


A very rough solution of equation (22) may be ean in the 
following way. 
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Put 
py ei == AS (23) 
Equation (22) now becomes 
Ag —" Oem Ae (24) 
The abscissa of the minimum of the left side is 
b 1/4 
m — hare ae 25 
x a( a ) (25) 
where a is a numerical coefficient. The ordinate of that minimum is 


Un = (08 — a) (> Ne (26) 


The point x, of intersection of the left side with the axis of abscissa 


1s: 
1b \v4 
=| — 24 
Lo ( a ) : (27) 


Very roughly we may approximate the portion yn of the curve 
(Figure 16) by a straight line passing through those points. The 
equation of this straight line is 


FIGURE 16 


b \14 b(a® — a) 
fo Oe 


A l-a 
Substituting that for the left side of (24) we find 


b \ 
r=-fa+(Z) ¢ (28) 
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For not too large masses, the second term of the right side of (28) 
prevails, and, remembering (13) we find: 

y, « M3, (29) 
From (29) and (20) we find 

Mz«M. (30) 


Relations (29) and (30) are fairly well satisfied by actual data, 
as seen from Figures 17 and 18. 


FicuRE 17. Frequency v, of breathing plotted against mass M of the animal 


on logarithmic scale. The straight line has the theoretically required slope of 
— 1/3. 1— rat; 2 — rabbit; 3 — cat; 4 — dog; 5 — sheep; 6 — bear; 7 — 
horse. Data from Crile (1941) and Oppenheimer and Pincussen (1926, p. 476). 


An interesting problem is presented in this connection by the 
lungs of birds. Their relative size is rather small, compared to those 
of mammals. But their structure and their mechanism of action are 
apparently more complicated (A. Bethe, 1925). The lung seems to 
be bathed by a continuous stream of air, due to a peculiar mechanism. 
From considerations of diffusion it can be readily seen that the ven- 
tilating efficiency of a bird’s lung is several times higher than that 
of a mammal. This permits a smaller volume, which is essential in a 
bird in order to reduce weight and: give room for the muscles operat- 
ing the wings. It may perhaps be possible in the future to show 
mathematically, by developing a more exact theory of the lungs, that 
a lung of the mammal’s type would be incompatible with flying and 
other functions of the birds, and that the structure of the bird’s lung 
is the simplest that satifies the necessary requirements. 

Somewhat similar considerations may be made about the heart. 


N. RASHEVSKY 53 


500 


200 


2 3 4 5678910 20 30 50 100 
—> M IN KG. 
FIGURE 18. Mass M, of lungs plotted against mass M of animal on loga- 


rithmic scale. The straight line has the theoretically required slope of 1. 1 — 
cat; 2 — cat; 3 — hare arctic; 4 — hare african; 5 — jackal; 6 — beaver; 
7 — coyote; 8 — dog; 9 — dog; 10 — greyhound; 11 — police collie; 12 — capi- 
bra; 13 — mountain lion; 14 — huskies; 15 — jaguar; 16 — leopard. Data from 
Crile (1941). 


Consider a capillary of length / and radius 7. Let that capillary 
supply with oxygen a region of radius R and length 1. Denoting by 
€ the average oxygen concentration in the capillary, assuming no ra- 
dial gradients, and denoting by c* the concentration of oxygen in 
the tissue, by h the permeability and by ag the rate of oxygen con- 
sumption per cm* per sec. proportional to the rate of energy metab- 


olism g , we have 


22h (GC — Cy — a Uh F* lag (31) 
or 
2rh 
Rk? =77 + —— (€— ¢"). (32) 
aq 


If é and c’* are considered as approximately constant, then for a 
given r equation (32) gives R. The higher q, the smaller RF, hence 
the larger the number of capillaries per cross section of tissue. Since 
q increases with decreasing mass, according to Rubner’s law, smaller 
animals should have a larger number of capillaries per unit cross sec- 
tion of tissue. This seems to be the case (Oppenheimer and Pincus- 


sen, 1925, p. 121). 
If v denotes the velocity of blood flow, then a7*vé is the amount 
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of oxygen entering the capillary at the arterial end. Considering for 
simplicity that at the venous end the concentration of oxygen is zero 
(the general case is treated similarly) we have: 


nr? vé = 2arlh(é — c*), (33) 
or, putting 
€—C—c. (34) 
UTC ZING, (35) 
or 
2lhe 
y= — (36) 
re 


From Poiseuille’s law we have for the pressure difference along 
the whole capillary: 


lv 
Ap =a—; 4=8u; uw — viscosity. (37) 
2 
Equations (36) and (37) give 
[? 
Ap=b—, (38) 
r3 
where 
2hac 
(— — (39) 
Cc 


Most of the blood pressure is across the capillaries. If we schem- 
atize the circulatory circuit, as shown on Figure 19; we see that (38) 


ARTERY 
an nee 
VEIN 


FIGURE 19 


gives approximately the average blood pressure. Approximately, in 
similarly built animals, the length | of the capillaries varies like the 
linear dimensions of the animal or like WM. If r is the same for all 
animals, then dp « M**. Actually Ap varies little with M , within a 
rather wide range of M’s. (Oppenheimer and Pincussen,. 1925, p. 
143). Hence, if Ap is approximately constant, we must have: 


pe met, (40) 


There may be different reasons for such a relation. They may be 
purely mechanical. A larger capillary may have to be thicker in order 
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not to break due to natural stresses, to which it is subjected by the 
tissues. 

One might think that an upper limit for r is set by the require- 
ment that the radial gradients should be very small. If ¢, and c? are 
the concentrations on the axis and at the periphery of the capillary, 
so that 


@. SP Gs 


é= ; 41 
5 (41) 
then we have, denoting by D the diffusion coefficient, 
2D CG ay Co) 
DTCs — Che (42) 
. 
or 
re h (43) 
CG, —C) es DS 
If 
Cz ro Co 
Os = GO 


should be below a small quantity «, we have 


This does not connect 7 with / in any way. 

To derive equations for the size and frequency of the heart, we 
may argue as follows: 

Assuming the beat volume of the heart to be roughly propor- 
tional to the actual maximum volume V,, and denoting by » the 


frequency, we have: 
Ve io Mq x iP", (44) 


according to Rubner’s law. 
The mass of the heart is proportional to V,, its linear dimensions 


to V./2. Hence the acceleration is proportional to V.1 ».2, and the 
inertial force to V.4/2 v2. The total force necessary to overcome the 
blood pressure p is proportional to pV,7/?. The mechanical force of 
the heart muscle is proportional to V.°”, their cross section. Hence de- 
noting by a, b and ¢ coefficients of proportionality, we have 


aV 2 v2 + bpV 22 =eV2", (45) 
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or , 
er 
V1 (46) 
a. 
Equation (44) solved for V., and introduced into (46) gives: 
As/4 
= (47) 
M2 
and 7 
M.= : (48) 
q At 


To what extent relations (47) and (48) are actually satisfied in na- 
ture is shown on Figures 20 and 21. 


— > Uc IN SECT! 
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Ficure 20. Heart frequency », plotted against mass M of the animal on 


legen seale. The straight ine. has the theoretically required slope — 1/3. 
1 — rat; 2 — rabbit; 3 — dog; 4 — sheep; 5 — tiger; 6 — lion; 7 — pig; 
8 — camel ; 9 — horse. Data from Crile (1941) and Oppenheimer and Pincus- 
sen (1925), p. 138). 


IX. Some Remarks on Unicellular Organisms 

If the principle of determination of structure by function is to 
hold generally, it should apply also to single cells. Physicochemical 
functions are more likely to be of importance here, than locomotion. 
The determining variables may be different chemical parameters. 
The rates of reactions are dependent on the concentrations of metab- 
olites, and the average concentrations of those are determined by the 
solutions of the diffusion equations, which involve the geometrical 
shape and structure of the system. The solution of the fundamental 
problem of the nucleus and cytoplasm may be sought in that direc- 
tion. If a cell of given mass is to be the seat of certain reactions, | 
which are quantitatively specified as to their rates, and if that cell is 
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FIGURE 21. Mass M, of heart plotted against the mass M of the animal on 


logarithmic scale. The straight line has the theoretically required slope 1. 
1 — cat; 2 — agouti; 3 — cat; 4 — hare arctic; 5 — hare african; 6 — jackal; 
7 — fox; 8 — rabbit; 9 — beaver; 10 — coyote; 11 — dog; 12 — capibra; 
13 — dog; 14 — greyhound; 15 — police collie; 16 — capibra; 17 — mountain 
lion; 18 — huskies; 19 — jaguar; 20 — leopard; 21 — hyena; 22 — lion; 
23 — lion. Data from Crile (1941). 


characterized by the presence of catalists, which control those reac- 
tions, then it may well happen that the necessary results may be ob- 
tained only if some catalists and some reactions are spatially segre- 
gated in a region which we call nucleus. A number of new problems 
in cellular biophysics are suggested by such considerations. 

The author is indebted to Professor Alston S. Householder for 
a critical reading and discussion of the manuscript of this paper. 

This investigation was aided in part by a grant from Dr. Wal- 
lace C. and Clara A. Abbott Memorial Fund of the University of 


Chicago. 
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SCOPE OF THE BULLETIN 

1. The Bulletin is devoted to publications of research in Mathe- 
matical Biophysics, as contributing to the physicomathematical foun- 
dations of biology in their most general scope. 

2. Papers published in the Bulletin cover physicomathematical 
theories as well as any other mathematical treatments of biological 
phenomena, with the exception of purely statistical studies. 

3. Mathematical studies in physics or in borderline fields in 
which a direct connection with biological problems is pointed out are 
also accepted. 

4, Emphasis is put upon the mathematical developments, but a 
description and discussion of experimental work falls also within the 
scope of the Bulletin provided that description or discussion is-made 
in close connection with mathematical developments contained in the 
same paper. 

5. Outside of the scope of the journal are papers of purely sta- 
tistical nature or papers concerned only with empirical equations. 


PREPARATION OF MANUSCRIPTS 

All manuscripts should be typewritten double space. Equations 
should be numbered on the right consecutively. Do not use such nota- 
tions as “equation 2a” or, “equation 5’” etc. References should be all 
given at the end of the paper, arranged in alphabetic order by authors 
and for each author chronologically, following exactly the style used in 
the Bulletin. In text, reference should be made by giving in paren- 
theses the name of the author followed by the year of publication. In 
case of several publications by the same author in the same year, use 
notations “1940a’’, ‘‘1940b’, etc. 

In writing equations, slanted lines should be used wherever pos- 
sible. 

Every paper is to be preceded by a short abstract. 

Drawings should be prepared in a professional manner on white 
paper or tracing cloth following as closely as possible the style used 
in the Bulletin. They should be drawn approximately twice the scale 
of the finished reproduction. Lettering should be made by means of a 
lettering guide. If desired by the authors, drawings may be made at 
the Editorial Office according to author’s instructions, at cost. 
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